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Abstract. We present a reanalysis of the He n Lyo- absorption towards the quasars HS 1700+6416 and HE 2347-4342 using new 
high S/N, optical observations. An alternative analysis method is applied, which fits the high quality, optical Hi data directly 
to the He n spectrum. The results are compared to those inferred from standard line profile analyses. This new method enables 
us to derive redshift scales characterizing the fluctuations of the column density ratio rj. We find n changing smoothly with 
redshift on typical scales of Az ~ 0.01 - 0.03 corresponding to 8 - 24 /r 1 Mpc comoving. The real length scales of variations 
of the column density ratio might be even larger, since part of the fluctuations may be caused by noise in the He n data and 
by systematic effects due to the applied method. However, variations on small scales cannot be ruled out completely. For n 
variations on scales of a few Mpc an amplitude of about ±1.5 dex cannot be excluded. The data shows an apparent correlation 
between low n regions and the presence of metal line absorbers, which corresponds to the more general correlation of low r\ and 
strong Hi absorption. Thermal line broadening is suggested as a probable explanation for this apparent correlation, since both 
fit methods would severely underestimate n for absorbers with A?hi <: 10 13 cnr 2 if the line width was dominated by thermal 
broadening. Indeed, lines located close to the cut-off of the b(N) distribution yield lower column density ratios compared to 
the whole sample, in particular if high density absorbers are considered. We argue that the apparent correlation of r\ with the 
strength of the H i absorption is caused by insufficient consideration of thermal broadened lines by the standard analysis. As 
unbiased value of the column density ratio, we find ^ ~ 80 in agreement with previous estimates. 
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1. Introduction 

The metagalactic ionizing radiation plays an important role for 
QSO absorption line studies probing the intergalactic medium 
(IGM). The IGM is highly photoionized and the UV back- 
ground governs its ionization state. Estimates of the intergalac- 
tic metallicity are highly dependent on the presumed ioniz- 
ing radiation (e.g. Agui rre et all2004T) . Nevertheless, the shape 
of the UV background is inaccessible to direct observation. 
Theoretical models based on the number, distribution, and 
spectral properties of underlying sources like AG N and star- 
burst galaxies have been computed bv e.g. | Ha ardt & Madau 
lll996ll200lh : lFardal et alJdl998l) : lBianchi et al.N200ll) . Efforts 
to derive the shape of the ionizing radiation from the ob- 
served column densities of different metal ions via photoio niza- 
tion c alculations have been made recently by Agafo nova et alJ 
J2005h . The observational studies suggest deviations from the 
shape predict ed by "classical" t heoretical models (see also 
iBoksenberg et all2003llFechner et al.|2006bl). Additionally, re- 
cent theoretical work JSchavd I20041 lMiralda-Escudel l2005) 
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pointed out that local radiation sources may be important for 
high density absorbers. 

Apart from column density ratios of metal line absorption, 
the Lya transition of He n (Aq = 303.7822 A) provides a further 
possibility to probe the UV background observationally. The 
column density ratio rj - Nh^i/Nhi depends on the photoion- 
ization rates of H i and He n and thus on the ionizing radiation 
at the H i ( 1 Ryd) and Hen (4Ryd) ionization threshold (e.g. 
iFardal et al.lll998|) . Furthermore, Hen probes the low density 
phase of the IGM, since He n is more difficult to ionize than H i 
due to lower fluxes and cross sections at its ionizing threshold 
and the faster recombination of He m compared to H n. 

Direct observations of Hen Lya absorption are very 
diffi cult due to the lim ited number of unabsorbed quasars 
dPicard & Jako bsen 119931) . Up to now, six objects have been 
observed, showing He n absorption in their UV spectra. Only 
two of them (HE 2347-4243 at z = 2.885 and HS 1700+6416 
at z — 2.72) are brig ht enough for high resolution spectroscopy 
with FUSE dKriss et . alJkOOll IShull et al ll2004| Iztieng et all 
120041 iReimers et al.ll2006UFechner et al.ll20Q6aK . The results 
from these data suggested that the UV background might 
fluctuate on small scales (Az ~ 0.001) with an amplitude 
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of several orders of magnitude (rj < 10 up to > 1000). 
Probable reasons for these fluctuations are the pro perties o f the 
OSOs like the spread in the spectral indices (e.g. lTelfer et alJ 
2002) an d finite lifetimes, lo cal density variations of the IGM 
jMiralda-Escude et alj gOOO) or filtering of radiation by radia- 
tive transfer effects jMaselli & Ferrarall2005l) . Furthermore, it 
has been shown that also the no ise level of the observed spectra 
contr ibutes to the scatter in rj jFechner et alJl20o"6al iLiu et alJ 
2006). 

Due to the limited quality of the FUSE data (S/N ~ 5) it is 
important to exami ne in detai l the ap plicability of the adopted 
analysis methods. Shull et al. (2004) applied an apparent opti- 
cal depth method, while all o ther studies are based on line pro- 
file fitting jKriss et alJ|200ll IZheng et~ai1l2004t lieimers et alJ 



Table 1. Observation diary of HE 2347-4342. All exposures 
use a 1" slit and 2x2 binning. The seeing conditions were 
very good with 0.5" - 1.0". 



2006). Fox et al. (2005) showed that the apparent optical depth 
method overestimates the column density in case of low signal- 
to-noise. Further more, it is only valid for unsaturated lines 
JSavage & Sembachll99lh. The applicability of line profile fit- 
ting was investigated by B olton et alJ 1120061) . who find that it 
leads to reliable results (but see discussion below), and the me- 
dian provides more confident results than the mean. According 
to their study, the fluctuations in rj are due to the small num- 
ber of quasars that contribute to the He n ionization rate at any 
given point. They also conclude that a small fraction of space 
have to be exposed to a harder radiation, since they failed to 
reproduce c olumn density ratios < 1, which are me asured by 
Zhen g et alJ J2004). A similar finding is reported bv lLiu et"afl 
(120061) . 

Here, we re-analyze the He n Lya forest towards HE 2347- 
4342 and HS 1700+6416 applying a new method, which fits the 
optical data directly to the FUSE spectra. This approach min- 
imizes effects of noise mimicking rj variations when applying 
the apparent optical depth method, and it avoids the subjectiv- 
ity of the line profile modeling, where the number of assumed 
components are due to arbitrary choice. In order to provide an 
alternative estimate of the scales on which rj fluctuates, we fol- 
low several approaches. Furthermore, we investigate the effect 
of line broadening on the inferred rj values. 

The paper is organized as follows: New high quality, optical 
data of the H i Lya forest towards HE 2347-4342 are presented 
in Sect. |2] In Sect. [3] we introduce the fit procedure (called 
spectrum fitting), which is applied to the data in Sect.|4] The re- 
sults, possible effects due to thermal line broadening, and their 
implications are discussed in Sect. [5] We conclude in Sect. [6] 

2. Data 

HE 2347-4342 was observed with UVES at the second VLT 
unit telescope (Kueyen) during 13 nights at two epochs in 2000 
and 2003. The details of the settings are presented in Tabled 
In 2000, the observations were carried out as a part of the Large 
Program "The Cosmic Evolution of the Intergalactic Medium" 
(1 16.A-0106A) in the UVES dichroic modes resulting in a to- 
tal wavelength coverage of ~ 3000 - 10400 A. Additional ob- 
servations were performed in 2003 in the blue with a central 
wavelength of 4370 A. The slit width of 1" and binning 2x2 
were chosen to match the setting of the Large Program data. 
The resolution is R « 45 000. 



date 


mode 


A cm (nm) 




2000/10/07 


DICHR#1 


346/580 


14400 


2000/10/08 


DICHR#1 


346/580 


7 200 


2000/10/08 


DICHR#2 


437/860 


7 200 


2000/10/09 


DICHR#2 


437/860 


7 200 


2000/11/19 


DICHR#2 


437/860 


7 200 


2003/06/05 


BLUE 


437 


7 200 


2003/06/06 


BLUE 


437 


10 800 


2003/06/07 


BLUE 


437 


6900 


2003/06/08 


BLUE 


437 


10 300 


2003/06/12 


BLUE 


437 


3130 


2003/06/17 


BLUE 


437 


6 600 


2003/06/24 


BLUE 


437 


6 600 


2003/07/06 


BLUE 


437 


3 700 


2003/08/05 


BLUE 


437 


3 244 



The data reduction was performed at Quality Control 
Garching using the U VES pipeline Data Reduction Software 
( Ball ester et all200ol) . The vacuum-barycentric corrected spec- 
tra were co-added. Thus, for the wavelength range 3800 - 
4950 A the total exposure time is 28.24 hrs resulting in a signal- 
to-noise ratio of S/N ~ 100 in the Lya forest. 

The UV data of HE 2347-4342 data have been taken with 
the FUSE s atellite in 2000. The data were first presented by 
Krissetal. (2001), and was newly reduce d and re-analyzed 
twice JShull et al.ll2004t IZheng et alJl2004 . In this study we 
use the spectrum of 1Zheng^tTlT ( l2004l) . and refer to their work 
for details of the data reduction. The resulting spectrum cov- 
ers the spectral range 900 - 1 188 A. The signal-to-noise ratio is 
S/N ~ 5onabin size of0.5Ainthe best part (1088 - 1 134 A) 
and slightly less elsewhere. The resolution is R « 20 000. 

A low resolution STIS spectrum was taken simultaneously 
with the FUSE observations in August and October 2000 using 
the gratings G140L and G230L. The spectra cover the wave- 
length range 1120 - 3 170 A and are used to extrapolate the 
continuum to the FUSE spectral range. We estimate the con- 
tinuum considering galactic extinction, Lyman limit absorp- 
tion, and the intrinsic spectral energy distribution of the QSO. 
The extinction curve is adopted from Carfejhetal. Jl989l) as- 
suming E(B - V) = 0.014 JSchlegel et alJ Il998b . In addi- 
tion, there is a Lyman limit system (LLS) on the line of sight 
atz * 2.735. From the Lyman break we estimate the op- 
tical depth r ~ 1.8 leadi ng to log/VLLS ~ 17.46 following 
Mcil ler & Jakobsenl Jl990l). which is in good agreement with 
the results of iRehiiersetai 1 dl997l) . However, the contributed 
optical depth of the LLS to the absorption in the FUSE spec- 
tral range is negligible (t < 0.08 at A < 1200 A) thanks to the 
v -3 dependence of the Lyman continuum opacity. Considering 
the LLS optical depth and the extinction law, we estimate the 
intrinsic QSO power law to have a spectral slope of a « -0.46. 

The UV spect rum of HS 1 700+64 1 6 taken with FUSE was 
first presented in iReimers et alJ d2006l). We u se the reduced 
and normalized data set from Fechn er et alJ (I2006al) . which 
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has a resolution of R as 20000 and a signal-to-noise ratio 
of ~ 5. The corresponding Hi Lya forest was observed with 
Keck/H IRES. For a de tailed presentation of the optical data we 
refer to lFechner et alJ J2006blal) . The resolution is R ta 38 000 
and S /W« 100. 

3. Method 

Previous studies of the Hen Lya forest are based on a 
profile-fitting proced ure dKriss et alJl200lt IZheng et alj r2004: 
iReimers et alJ 120061) or the apparent optical depth method 
fchull et alJl2004 . Both methods suffer from different short- 
comings. The apparent optical depth m ethod leads to reli- 
able results only if lines are optically thin llSavage & Sem bach 
Il99lh . In case of the Hen forest many profiles are saturated. 
Thus, the validity of the apparent optical depth method is lim- 
ited. Furthermore, the results strongly depend on the noise, 
which may mimic fluctuations of r\. 

In case of fitting line profiles to the data the H i spectrum is 
modeled by Doppler profiles. The derived line parameters are 
used to fit the He n data as well. In doing so, the line widths 
have been assumed to be dominated by turbulent broadening, 
i.e. £>Hen = t>Hi- This assumption is also made implicitly in the 
apparent optical depth method, and is inherent in our spectrum 
fitting method as well. However, fitting line profiles need the 
subjective decision how many components are used to model 
the observed features. Particularly, this problem gets significant 
for weak lines and complicated blends, where components may 
get lost or the decomposition is insufficient. Problems in the H i 
model may cause a poor description of the corresponding He n 
absorption, where the features are stronger but the data quality 
is worse. Therefore, apparent variations in 77 may be artifacts of 
an inadequate decomposition. 

With the purpose to avoid a strong dependence on the noise 
as well as the subjectivity of a model, we develop a new method 
to estimate 77. It is based on the idea to compare the He n and 
the H i Lya forest directly without the need to provide a model 
of the data and without using a pixel by pixel (or bin by bin) 
technique which would be affected by the noise level of the 
data. Therefore, we fit the whole H i spectrum to the H e n data 
by scaling it by 77 = 4 • r He n /r H i < Miralda-Escudel 1 993h. which 
is then the only free parameter. 

Since an unbiased H 1 spectrum is needed, we identify the 
metal line systems in the optical data, fit Doppler profiles to the 
spectral range of interest and subtract them from the spectrum. 
Furthermore, we confine ourselves to the pure Lya forest, and 
only the data between Lya and Ly/3 emission are suitable. In 
case of HE 2347-4342 these are wavelengths > 4000 A in the 
optical and > 1000 A in the UV. Strong Ovi absorption aris- 
ing from a complex system associated with the QSO is present 
shortward of ~ 4050 A. Due to the complexity of the system 
and severe blendin g with forest lines (for a more detailed dis- 
cussion see iFechneret alJ 2004i. we restrict the redshift range 
to 2.33 < z < 2.91. 

In preparation of a fit, the wavelength scales of the 
datasets are aligned scaling the optical data by a factor of 
1215.6701 A/303.7822 A * 4.00178. Furthermore, the resolu- 
tion of the optical spectrum is changed to match the lower reso- 



lution of the FUSE data by a convolution with a Gaussian. The 
method assumes implicitly that line broadening is dominated 
by turbulent motion and the line widths are the same for H 1 and 
He n. Since the operation of convolution is commutative, the 
width of the Gaussian has to be <x = (c^ =20000 - (^^qqq) 1 ^ 2 - 
In order to perform the actual fit, a wavelength portion is 
chosen and log 77 is estimated to minimize 



, 1 (/He II, i ~ exp (- 1 • T H ,, i)) 2 

x = - ) .- 



°Hbn.i + <7 Hi.j 



(1) 



where fii eu and thi = -In /m are the normalized observed 
flux and the optical depth of Hen and Hi, respectively, with 
the uncertainties <x He n and <x H : ; « is the number of pixels in the 
chosen wavelength interval. The 1 cr uncertainty of the fitted 
log 77 is estimated finding the values at which y 1 — X 2 ■ ± 4f 2 
with Ay 2 = 1 .0. We restrict the parameter space to values 0.0 < 
log 77 < 4.5. 

The idea of fitting the spectra directly is illustrated in Fig. 
□ for the spectral range 1088.0 - 1 133.5 A of HE 2347-4342, 
where the data quality of the UV spectrum is best. Each panel 
shows the normalized FUSE data in comparison to the aligned, 
metal line corrected Hi spectrum scaled by log 77 = 0.0 (i.e. 
unsealed), 1.0, 2.0, and 3.0, respectively. Obviously, some re- 
gions are consistent with log 77 as 1 (e.g. 1112 - 11 17 A) or 
log 77 « 2 (e.g. 1096 - 1103 A). Whereas scaling the Hi spec- 
trum by log 77 3 leads to almost complete absorption with 
few transmission windows, which is apparently not consistent 
with the Hen observations in the presented wavelength range. 
The overall fit of this spectral range leads to log 77 = 1.81^ ^ 
with^- 2 = 3.34. However, the considered portion corresponds 
to Az ~ 0.15, and obviously, 77 varies on smaller scales. 

The definition of a scale in redshift, on which 77 is sup- 
posed to appear constant, is the most crucial point. The required 
length scale is given by the number of pixels n that are consid- 
ered. If n is chosen too small, the result is severely affected by 
the noise level of the FUSE data. If n gets very small, the pro- 
cedure approaches the apparent optical depth method. On the 
other hand, possible small-scale fluctuations will be smoothed 
and therefore remain undetected if the length scale is chosen 
too broad. 

In order to investigate possible effects on the results, we 
follow four different approaches, (i) constant binning, (ii) bin- 
ning "per eye", (iii) binning by ^-estimation, and (iv) aver- 
aging around each pixel. We introduce these procedures in the 
following. 

(i) Most simply, bins of a fixed size are adopted. In this 
case, the bin size (we use 0.5 and 1.0 A) and the positioning of 
the bins are free parameters. Therefore, fluctuations are inves- 
tigated only on the given scale. 

(ii) In order to find a length scale, which is not fixed a pri- 
ori, we define the fit intervals by eye. Intervals are fitted suc- 
cessively starting at the low wavelength edge of the data. The 
size of the intervals is chosen to be as large as possible to show 
no obvious fluctuations of 77. The advantage is that apparent 
changes of 77 as illustrated by Fig.^can be considered. On the 
other hand, the choice of the size of the intervals is made very 
subjectively. 
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lQg7)= 0.0 




Fig. 1. Comparison of the observed He n data towards HE 2347-4342 to the metal subtracted, convolved H i spectrum scaled by 
log 77 = 0.0, 1.0, 2.0, and 3.0. Apparently, there are regions modeled well by log?/ » 1.0 (e.g. 1112 - 1117 A) or log 77 » 2.0 
(e.g. 1096 - 1103 A and ~ 1125 A), respectively, while log 77 « 3.0 would lead only to sparse regions with non-zero UV flux, 
which is inconsistent with the observations in this redshift range. Features that appear nearly overestimated by log 77 = 0.0 (at 
~ 1 104 A, 1 108.5 A, and 1 121.5 A) correspond to weak metal line systems (see text). 



(iii) The more objective way is to start with a certain num- 
ber of pixels at the low wavelength limit of the spectrum and 
estimate the x 1 of the fit. More pixels are added to the fit in- 
terval as long as the^- 2 of the fit decreases. If the^- 2 increases 
again, the interval is defined and the next evaluation starts. Test 
calculations show that in the case of HE 2347-4342 at least 15 
pixels, corresponding to 0.375 A have to be used for the first 
evaluation to avoid confusion due to the noise. This means the 
smallest detectable scales have a size of Az ~ 0.001. 

(iv) Alternatively, we average over an interval of a certain 
size centered on every pixel. This procedure corresponds to a 
pixel by pixel estimation smoothed by averaging over a certain 
number of adjacent pixels. However, the results again depend 
on the number of pixels used, i.e. the degree of smoothing. A 
good choice appears to be 40 pixels corresponding to 1 .0 A. 

Another crucial point is the continuum estimation of the 
optical data. If the continuum is located too low, absorption 
of H 1 will be removed and the fitted 77 value will be overesti- 
mated, while 77 will be underestimated if the continuum level 
is too high. In order to determine an optimal continuum we 
perform the continuum fit simultaneously with a line profile 
fit u sing the line fitting program CANDALF (also used e.g. 
bv tFechner et alJ l2006a). The fit program models the contin- 
uum by Legendre polynomials. The simultaneous estimation 
of continuum and line parameters is supposed to lead to a more 
accurate continuum level than an a priori normalization. Using 



the formal errors of the Legendre parameters a statistical un- 
certainty of the continuum can be estimated. We find a statis- 
tical error of less than ~ 1 %. Additional systematic uncertain- 
ties due to the considered number of components are hard to 
quantify. Because of the exceptionally high signal-to-noise ra- 
tio of the optical data a realistic estimate of the systematic error 
should be < 1 %. 

Test calculations show that a continuum error of 1 % leads 
to an uncertainty in 77 of +0.1 dex. This value is small in com- 
parison to the fit errors due to the noisy He n data which are on 
average +0.4 dex or even larger (see below). Therefore, possi- 
ble uncertainties introduced by the continuum fit will be ne- 
glected in the following. However, the presented errors of log 77 
might be about 3 % larger due to the uncertainties in the con- 
tinuum estimation. A more conservative estimate of the contin- 
uum fit error (~ 2 %) would result in an additional uncertainty 
in 77 of roughly 0.2 dex which would enlarge the given errors by 
about 12%. 

4. Results 

Applying the four scale-defining methods described above, the 
data are analyzed. For comparison, we also perform a Doppler 
profile analysis using the line fitting program CANDALF. The 
resulting distributions of 77 with redshift in case of HE 2347- 
4342 are presented in Fig. |2] The panels show the results for 
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Fig. 2. Distribution of rj with redshift towards HE 2347-4342 defining the fit interval by eye (upper panel) or using a^ 2 procedure 
(middle panel) as well as the result from profile fitting (lower panel). For a better orientation, the histogram-like line in each 
panel indicates r\ estimated with a constant 1.0 A bin size. The redshifts of metal line absorption systems are marked as arrows 
in the upper panel. 



binning by eye, binning using the x 1 method, and from the 
Doppler profile analysis. For a better orientation, the result for 
constant binning of 1 .0 A bin size are shown as overlay in each 
panel. 

The low S IN of the FUSE data leads to large error bars, 
regardless of how the fit interval is defined. At longer wave- 
lengths (> 1135 A corresponding to z > 2.73), instead of 
a forest structure, the Hen absorption gets "patchy". This 
means that regions of high opacity and nearly complete ab- 
sorption change with opacity gaps and continuum windows 
jReimers et alJll997t ISmette et ali 2002). Unlike the Lya for- 
est, the absorption in the patchy zon e is thought to aris e 
from a medium not yet fully reionized jReimers et alJ fT997). 
Therefore, different results may be expected and we distin- 
guish between the Lya forest visible in Hen and the patchy 
absorption towards HE 2347-4342. The slightly less redshifted 
QSO HS 1700+6416 probes only material that is already com- 
pletely reionized. This means, solely Lya forest absorption is 
detected. In the following, we first present the results from the 
Lya forest in the spectrum of HE 2347-4342 and thereafter 
HS 1700+6416. 



4.1. The Lya forest towards HE 2347-4342 

As Lya forest towards HE 2347-4342 we consider the wave- 
length range < 1134 A or z < 2.735, respectively. As can be 
seen from Fig.|2] the rj values in the redshift range 2.58 < z < 
2.73 are characterized by comparably small error bars. A close 
up to this range, where the data quality of the FUSE spectrum 
is best (1088 - 1134 A), is shown in Fig. [3] In this figure the 
results from all different methods concerning the choice of the 
scale are presented. The fact that we find similar distributions 
of rj with redshift for all scale definitions (panels a - d), gives 
confidence that the choice of the length scale has a minor in- 
fluence on the results and can, in principle, be derived from 
the data. Even the column density ratios of the independent 
profile-fitting analysis are distributed similarly (Fig. [3^). This 
is noticeable in particular in the redshift regions 2.60 - 2.62, 
2.645 - 2.66, and 2.68 - 2.70. 

Extremely large error bars of the column density ratios ob- 
tained by profile fitting appear to be correlated to extreme 77 
values or large uncertainties in the corresponding results of the 
spectrum fitting procedure. Obviously, there are problems with 
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Fig. 3. Redshift distribution of 77 for the best part of the FUSE 
data (1088 - 1 134 A) of HE 2347-4342 using different fit inter- 
val definitions, (a): constant binning with bin sizes of 0.5 A 
(thin line) and 1.0 A (thick line; corresponding to the refer- 
ence distribution in each panel of Fig.|2}. The uncertainties are 
omitted for a clearer presentation, (b): fit interval definition by 
eye. (c): bin size estimated by ax 2 procedure, (d): estimation 
of rj per pixel averaged over 1.0 A, i.e. 40 pixels. The dotted 
line represents the 1 cr uncertainties, (e): result from profile fit- 
ting. The dashed lines in each panel indicate the 77 distribution 
yielded with the x 2 procedure (panel c) for better comparabil- 
ity. All presented errors specify 1 <x uncertainties. The positions 
of metal absorption systems are indicated as arrows in panel a. 

the data at ~ 1093 A, where several noise peaks in the FUSE 
data mimic an absorption structure different from that of Hi 
(see also Fig.^. Strong saturation, i.e. nearly complete absorp- 
tion over roughly 4 A, lead to large error bars at ~ 1 127 A. A 
similar configuration produces the large uncertainties of the re- 
sults from line profile fitting at ~ 1 105 A although the spectrum 
fitting results remain nearly unaffected. 

The median value for this best region is log 1.75 in case 
of the spectrum fitting method, similar for all adopted scales, 
and 1.80 in case of the profile fits considering also lower lim- 
its. TTrjsj£_sHghtlyless than the mean value of 77 — 80 found 
bv lKriss et alJ J200ll) . The 1 cr scatter related to the average is 
roughly 0.7 dex and slightly larger (~ 0.8 dex) for the whole 
Lya forest observed. Furthermore, almost all derived 77 values 
are below log 77 = 3.0. Since very high 77 values are produced 



Table 2. Metal line systems with z > 2.3 in the spectrum of 
HE 2347-4342. 



z observed ions 



2.3132 a 


Cm, 


Civ, Sim, Siiv 


2.3475 


Cm c 


, C iv, O vi b 


2.3741 


Ovi b 




ZAioZ 


f TTT C 
L, III 


, U IV, U VI 


2.4887 


Civ, 


Ovf 


2.6346 


Civ, 


Ovi 


2.6497 


Cnf 


, C iv, O vi 


2.6910 


Ovi a 




2.7121 


Civ, 


Ovf 1 


2.735 d 


Cii, 


Cm, Ovi, Sim, Siiv 
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a apparently related to a QSO close to the line of sight at z = 2.31 

(Jakobsen, private communication) 
b only one doublet component detected 
c uncertain detection due to blends 
d Lyman limit system 
c complex associated system 



by soft radiation such as stellar sources, this might indicate that 
even though galaxies contribute to the ionizing radiation, they 
are not dominating the intergalactic UV background at the ob- 
served epoch. 

Values of 77 significantly below the mean are visible at He n 
wavelengths around 1090 A, 1104 A, 1108 A, and 1121 A in 
Fig-El in addition a local minimum is located at 1 128 A. Four 
of these minima coincide with metal line absorption systems 
at z = 2.6346, 2.6497, 2.6910, and 2.7121 (marked with ar- 
rows in Fig. [3^ and summarized in Table |3 associated with 
saturated Hi absorption features. The 77 minimum at 1090 A 
is due to a Hi line (logN ~ 13.0), which seems to mismatch 
the absorption features in He 11 (see Fig. probably an ef- 
fect of the poor data quality at this spectral region. The metal 
line systems show weak C iv absorption (log N < 12.8), promi- 
nent Ovi features, and no low ionization species. Since the 
O vi features are located in the Lya forest, their identification 
might be questionable due to blending with forest lines in at 
least one doublet component. Since the systemic redshift can 
be fixed by the unblended C iv, the existence of those systems 
is confirmed. However, the system at z = 2.6910 shows no 
Civ absorption and the /11038 component of Ovi is located 
in the Ly/3 absorption trough of the LLS at z ~ 2.735. Thus 
the only evidence for metal line absorption is provided by the 
A 1032 component of O vi. Since the detected feature is narrow, 
b = (11.06 + 0.81)kms _1 , it certainly originates from a metal 
ion although an alternative identification cannot be ruled out 
completely. 

Further apparent minima of 77 are found at lower redshift. 
The position of metal absorption systems are indicated as ar- 
rows in the upper panel of Fig. [3 The four metal line systems 
with 2.33 < z < 2.50 showing only highly ionized species coin- 
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Fig. 4. Distribution of 77 in voids and filaments derived from 
binning with constant bin size (1 A; left panel) and line profile 
analysis (right panel) in the Lya forest region (z £ 2.73) to- 
wards HE 2347-4342. Following Ishull et alJ fc004 . voids are 
defined as bins with a mean opacity in hydrogen of thi < 0.05. 
Concerning the results from the line profile-fitting procedure, 
we consider lines with log Nri > 13.0 to arise from filaments. 

cide again with minima of 77. The system at z = 2.4887 appears 
to be located in between of two low 77 regions according to the 
upper panel of Fig. |2] However, referring to the lower panels, 
only the minimum at z ~ 2.48 may be real. 

Additional low 77 regions are visible at z ~ 2.39 and 2.53. 
The latter might be affected by the extremely high noise in the 
corresponding He 11 wavelength range, while the first one is an 
apparent void in H 1. Within Az a 0.063 (2.368 < z < 2.430) we 
only detect lines with logA^i 5; 13.5. At the central position 
of the void (z « 2.398), two stronger lines (logA^i ~ 13.5) 
arise. We suspect that this void could be due to a quasar close 
to the line of sight, whose hard radiation field ionizes H 1 and 
Hen resulting in the observed low 77. Fitting the range 2.388 < 
z < 2.401 yields log 77 = 1.72*} ^. The appearance of the fit 
suggests fluctuations on scales smaller than Az « 0.013 also 
reflec ted in the large 1 <x errors. 

Shull et al1lEo 04) found a correlation between the strength 
of the Hi absorption and 77, in the sense that higher 77 val- 
ues are measured in Hi voids. This finding was confirmed 
by iReimers et ail feOOd) towards HS 1700+6416. We re- 
investigate this point estimating the mean Hi opacity of the 
considered fit intervals. The distribution of 77 for constant bins 
of 1.0A size of high (t H i > 0.05) and low (r H i < 0.05) 
opacity is shown in the left panel of Fig. 0] Again, only the 
Lya forest range is considered. The mean values are log 77 = 
1.53 + 0.79 and 1.98 + 0.71 in the high and low opacity sam- 
ple, respectively. For comparison, the corresponding distribu- 
tion obtained from the profile-fitting procedure is presented in 
the right panel of Fig.0] The optical depth of t = 0.05 at the 
central point of the Doppler profile for a H 1 Lya line of the 
width b = 27.0 km s~' is reached at \ogN = 12.25. In Lya 
forest statistics, a void is defined as a region without any ab- 
sorption features st ronger than log N h , ~ 13.5 over a certain 
comoving size ("e.g. lKim et aljEoOll) . As a compromise and in 
order to maintain a sufficient number of high density absorbers, 
we consider lines with log Af < 13.0 as low density absorbers. 
Furthermore, lines with column densities < 10 13 cirT 2 seem 
to remain unaffected by effects of thermal line broadening (see 



Sect.[5]for more details). For the low density sample we find the 
mean value log 77 = 1.90 + 0.74. Absorbers with logN > 13.0 
lead to log 77 = 1.21 + 1.50. Considering different samples of 
strong H 1 absorbers, we find that the mean value of log 77 de- 
creases with increasing line strength, either defined by optical 
depth or H 1 column den sity. This corresponds to the sugges- 
tion of IShull et alJ J2004I) that low density absorbers correlate 
with high 77 values. It is also consistent with the apparent cor- 
relation of metal line systems with strong Hi absorption and 
low 77 values. We suspect that this correlation is, at least partly, 
due to thermal line broadening, which lead to an underestima- 
tion of the Hen column density for absorption features with 
logA^Hi ^ 13.0 if turbulent line widths are assumed. We will 
return to this point for further discussion in Sect. [5] 

4.2. HS 1700+6416 

For the analysis of the He n data of the quasar HS 1700+6416, 
we follow the same method as described in Sect. [3] 
Additionally, t he metal lines predic ted to arise in the FUSE 
spectral range ( Fe chner et all2006bl) are considered during the 
optimization procedure modifying the ^-estimation (Eq.Q: 

2 1 (/He II, i - exp (- \ • T H 1, i ~ T met> ,)) 



This strategy is used, since due to the noise the results are more 
stable if the UV metal lines are considered as part of the fit 
instead of removing their optical depth from the FUSE data. 

Fig. [5] presents the redshift distribution of 77 towards 
HS 1700+6416 derived by spectrum fitting defining the fit 
interval by eye (upper panel), using a y 2 procedure (middle 
panel), or adopting profile fits (lower panel). Metal lines and in- 
terstellar absorption hav e been considered whe n fitting Doppler 
profiles (for details see iFechner et alJl2006ah . In the redshift 
range 2.60 < z < 2.75 there is apparently only little scatter in 
77. Indeed, the whole range (2.5815 < z < 2.7527, correspond- 
ing to Az = 0.1712) can be fitted with log 77 = 1.98+° f v which 
is in good agreement with the mean and median values result- 
ing from the applied methods. The low value of 77 at z w 2.735 
is due to continuum in H 1 and He n. In this case, the applica- 
bility of the spectrum fitting method is limited and it tends to 
produce low 77 values. 

Stronger fluctuations are visible in the lower redshift range 
although the search for a redshift scale using ^-estimation 
leads to smooth variation around log 77 ~ 1.8 at z = 2.479 + 
0.066. Fitting the whole range yields log 77 = 1.78!°;^. High 
77 values are apparently found at z « 2.41. This peak is most 
likely due to strong interstellar Cn A 1036 absorption, which 
is not considered in the spectrum fit procedure. Furthermore, 
strong Civ absorption from the Lyman limit system at z = 
2.3155 is expected in the UV, which is probably underpre- 
dicted (see Fechner et al. 2006b). Therefore, we overestimate 
the amount of absorption due to He 11 leading to unrealistically 
high values of 77. Considering the interstellar Cn absorption 
should reduce the derived values. Compared to the result from 
the Doppler profile fits, which includes the galactic C 11 line, 
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Fig. 5. Redshift distribution of 77 towards HS 1700+6416 according to the spectrum fitting method defining the fit interval by 
eye (upper panel) or using a x 1 procedure (middle panel) as well as the result from profile fitting (lower panel). For a better 
orientation, the histogram-like line in each panel indicates 77 estimated with a constant 1.0 A bin size. The redshifts of metal line 
ab sorption systems ar e marked as arrows in the upper panel, the position of galaxies close to the line of sight from the compilation 
of lShaplev et alJ d2005l) are indicated by thickmarks as well. 



no apparent high 77 peak is visible in agreement with the above 
arguments. 

Close to the metal line system at z — 2.4965, another re- 
gion of increased 77 can be seen, where the redshift distribution 
yielded by profile fitting resembles that found with the spec- 
trum fit. The most prominent Hi Lya feature is the central 
component in a region (2.485 < z < 2.504) of weak Hi ab- 
sorption, i.e. a void. In the corresponding He 11 spectral range 
strong absorption is detected. However, several data points of 
very low 77 in the range z ~ 2.50 - 2.5 1 derived from the profile 
fits reflect a major problem in the interpretation of the results 
concerning the zero flux level of the data. The H 1 features at 
this position appear stronger than the corresponding He 11 fea- 
tures. The same effect is also visible at z ~ 2.31, 2.435, and 
2.46, where the Hi absorption at the two first redshifts is par- 
ticularly strong, since they are the Lya features of Lyman limit 
systems. As can be seen from Fig. |6] the strong He 11 absorp- 
tion features at z ~ 2.435 appear to be saturated, even though 
the line cores do not reach zero. Therefore, conclusions drawn 
from the low redshift portion of the FUSE data are only pre- 



liminary. However, in the longer wavelength part of the data 
(> 1087 A, above the FUSE detector gap), no indications for 
an incorrect zero flux level are found. Thus, results from the 
range z > 2.58 can be assumed to be more reliable. 

In the spectrum of HS 1700+6416 25 meta l line systems 
can be identified (e.g. IVogel&Reimerslll995l iKohler et alJ 
Il996l iFechner et al.ll2006bl) . 13 of them are at redshifts z > 



2.29, i.e. their He 11 absorption is located in the observed FUSE 
spectral range. However, three are affected by the Hi Ly/3 air- 
glow emission or the detector gap. From the remaining sys- 
tems, only three arise in the reliable longer wavelength re- 
gion. All these are associated to the QSO showing absorp- 

J 1 1 n "—i 

tion features of O vi and N v ( Simcoe et al. 2002; Fechner et al. 

2006b). The system at z — 2.7443 is blueshifted with respect to 
the quasars' emission redshift a nd shows multicomponent O vi 
absorption dSimcoe et alJl2002l) . Since these systems are cer- 
tainly affected by the quasars' radiation, they are highly suit- 
able to investigate the effect of hard ionizing radiation on 77. 
Nevertheless, no indications for a proximity effect are notice- 
able in the distribution of 77 in Fig. [5] A slight dip might be 
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Fig. 6. Absorption features of H i Lya (thick line) and He n Lya 
(thin line) at z ~ 2.435. The profile of the Hen features indi- 
cates saturation, even though the flux does not reach zero. The 
dotted line represents the error of the He n flux. The narrow line 
in between both the saturated H i troughs is Si n of the LLS at 
z = 2.3155. 



present at z ~ 2.71 close to the metal line systems also visible 
in the profile-fitting result. Considering the blueshifted system, 
we do not find any hints for a depression of 77. 

The absence of a proximity effect in 77 seems to be in con- 
flict with the presence of associated metal line absorption obvi- 
ously exposed to hard quasar radiation. However, the observed 
redshifts of the metal line systems may not be solely cosmo- 
logical. This means, the absorbers may be located very close 
to the QSO (~ lOkpc), ejected from the central region with 
high velocities. The system redshifted with respect to the QSO 
emission z em ~ 2.72 is supposed to trace infalling material, i.e. 
the distance to the central engine is unkno wn. Due to the short 
duty cycle of quasars (~ 10 6 - 10 7 yr, e.g. lako bsen et al J20031 
Schirb er et al . 2004), the assumption of photoionization equi- 
librium might not be fulfilled for these systems. In ad dition, 
HS 1 700+6416 is variable in the intrinsic EUV jReimers et alJ 
120051) . However, short time variations are not supposed to cause 
deviations from ionization equilibrium. 

However, there seems to be an apparent correlation of metal 
line systems and regions of low 77 at redshifts < 2.55 (see Fig. 
|5J, which have to be interpreted carefully due to the discussed 
un certainties in th e zero fl ux level. 

iFechner et all (|2006b) found that the majority of the 
metal line systems with z > 2 towards HS 1700+6416 can 
be modelled con vincingly using a modified version of the 
lHaardt & Madaul J200ll) UV background. The modification of 
the shape of the ionizing spectrum is a shift of the break usually 
located at 4Ryd to lower energies (3 Ryd, HM3). The column 
density ratio 77 expected for this radiation field is 150 - 190, de- 
pending on the exact position of the break. If this shape of the 
ionizing radiation is representative for the IGM in general, the 
expected column density ratio would be log 77 ~ 2.2, consistent 
with the value found by fitting the H 1 spectrum to the upper 
FUSE spectral range. The metal line systems, however, probe 
the high density IGM. Therefore, according to the analysis of 
metal line systems, at least the high density absorbers are ex- 
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Fig. 7. Dependence of 77 on the strength of the H 1 absorption as 
derived from artificial absorption features in case of the spec- 
trum fitting method (upper panel) and the profile-fitting pro- 
cedure (lower panel). For several column densities of Hi and 
bui = 27.0 km s~', line profiles of Hi and Hen are calculated 
with 77 = 80 and bu eu = Z?hi (turbulent broadening, filled cir- 
cles) or />Hen = 0.5 ■ bii! (thermal broadening, open circles), 
respectively. The artificial data was analyzed assuming pure 
turbulent broadening. The dashed line indicates the presumed 
value of 77. 

pected to show high 77 values, which contradicts the results ob- 
tained by measuring H 1 and He 11 directly. We suspect that this 
conflict can be solved if thermal line broadening is taken into 
account. The distribution of the Doppler parameters derived 
from line profile fitting has its maximum at Z?hi ~ 23kms~ 1 . 
If this line width is interpreted to be completely thermal, the 
corresponding temperature is roughly T « 3 • 10 4 K, which i s 
close to the valu e 2 • 10 4 K estimated by Rico tti et all yOOO). 
lLai et al.l J2006I) even predict temperature fluctuations due to 
inhomogeneous He 11 reionization with temperatures as high as 
3 ■ 10 4 K. Thus, contributions from thermal broadening to the 
line width are expected to be important at least for part of the 
absorbers. 



5. Discussion 

The analysis of the He 11 data of HE 2347-4342 revealed an ap- 
parent correlation between regions of low 77 and the presence 
of metal line absorption systems. This finding is confirmed to- 
wards HS 1700+6416, albeit systematical uncertainties possi- 
bly affect the results of this line of sight. In order to investigate 
biases caused by the assumption made for the investigation, we 
analyze artifical data. Since turbulent line broadening is implic- 
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itly assumed in the spectrum fitting procedure, we consider two 
sets of profiles dominated by thermal and turbulent broadening, 
respectively. 

The artificial lines are computed as Doppler profiles with 
z = 2.6194, £>hi = 27.0 km s~' and different values of the 
Hi column density with 11.5 < logA^i < 15.5. The Hen 
column density is logA^Hen = log 77 + logA^i with 77 = 80. 
Two He 11 features are computed for each column density rep- 
resenting both the extremes of line broadening, pure turbulence 
'-'Hen = bni and thermal broadening buen = jbni- The gener- 
ated line profiles are convolved to match the resolution of the 
observed optical (R * 45 000) or UV spectra (R * 20000), re- 
spectively. Noise is added to have S /N ~ 100 in case of H 1 and 
~ 5 in case of He n. 

The simulated lines are analyzed twice, fitting the spectrum 
and Doppler profiles. Fig. shows the inferred 77 as a function 
of the underlying H 1 column density. Both fitting procedures 
severely underestimate 77 for strong Hi lines if the assump- 
tion of turbulent broadening is incorrect. Here we consider the 
extreme of pure thermal broadening, where the Doppler pa- 
rameter of He 11 is half the value of H 1. In reality, the re will 
be a m ixture of thermal and turbulent broadening. Z heng etaD 
(2004) even find that turbulent broadening is dominant. This 
finding has been confirmed bv lLiu et"aD (2006) who compared 
the o bservations to cosm ological hydrodynamical simulations. 
Also lBolton etail £2006) argue that thermal broadening is sup- 
posed to be negligible, since Hubble broadening is expected to 
be the dominant contribution to the width of Lya forest lines. 
However, the stronger the contribution of thermal broadening, 
the more the inferred 77 will underestimate the actual value. 

If turbulent broadening is indeed dominating, line profile 
fitting recovers well the underlying 77. In this case, the artificial 
lines are computed as Doppler profiles. Therefore, profile fits 
are supposed to recover the correct values, in particular since 
only a single line is considered. A more detailed exploration 
would require hydro dynamical simulat ions. Such an investiga- 
tion is presented by iBolton et al.l Il2006h exploring the appli- 
cability of line profile fitting. They find that 77 is infe rred accu- 
rately using this technique (but see discussion below). ILiu et alJ 
( 2006) additionally investigate the effect of thermal broadening 
and peculiar velocity fields using hydrodynamical simulations. 
They conclude that the absorption features are predominately 
turbulent broadened but fluctuations of the temperature and the 
velocity fields contribute to the scatter in 77. 

According to our synthetic line profiles, the spectrum fitting 
method exhibits deviations from the underlying 77 value even in 
the case of turbulent broadening. The reason is that we use the 
relation THen = 77/4 • thi to scale the H 1 spectrum. If the H 1 ab- 
sorption features are saturated, the accuracy of the optical depth 
measurement is limited. Though the correct 77 is recovered 
within the 1 <x error bars, the root-mean-square deviation is 
roughly 0.16 for column densities log Nhi > 13.5. Considering 
a sample of synthetic line profiles with log 77 = 1 .0, the spec- 
trum fit method recovers the underlying column density ratio 
more accurately in the case of stronger H 1 features. 

With this information about potential methodical problems, 
there are three possible interpretations of the apparent correla- 
tion between low 77 and strong Hi absorption: (1) a methodical 



bias, (2) an effect of thermal line broadening, or (3) physical 
reality. In the following, we will discuss these interpretations 
and their implications. 

5.1. Methodical uncertainties 

Due to methodical effects 77 might be systematically underesti- 
mated for strong Hi absorbers, since tHeii = 77/4 • thi leads to 
accurate results only for unsaturated lines. Since metal lines 
arise particularly in strong Lya systems, this would lead to 
the apparent correlation. At first glance only the results from 
spectrum fitting should be affected i n this way. In accord ance 
with Fig. 0and the investigation of IBolton et alJ J2006I) pro- 
file fits recover the underlying 77 values. Howeve r, the scatter 
plots o f the column density ratio 77 presented by B olton etafl 
(2006, their Fig. 4) show a slight tendency to smaller 77 values 
in the case of higher H 1 column densities, even for a uniform 
UV background. Due to the power law behaviour of th e col- 
umn density distribution function (e.g. Kirkman & Tytle rl 19971 
iKim et aDl200ll) . high column density absorbers are less nu- 
merous than weaker lines. Thus, strong absorbers represent 
only a minor fraction of the investigated line sample with a 
negligible contribution to the overall statistics. Therefore, it is 
necessary to extend the analysis of simulated data with the pur- 
pose to obtain a sufficiently large statistical sample in order to 
investigate high column density absorbers (logA^i k, 14.0) in 
more detail. 

5.2. Thermal broadening 

Assuming that the profile-fitting procedure recovers the col- 
umn density ratio accurately, the apparent correlation between 
metal line absorption and low 77 should be no methodical arti- 
fact. The reason is that the correlation is also recognizable in 
the profile-fitting results (Fig . [3 and [5j . Then, Fig. suggests 
that the line widths of Hi and He 11 features with associated 
metal line absorption might be dominated by thermal broad- 
ening. In general, lines with Nhi S 10 13 cm~ 2 if dominated 
by thermal broadening, would produce extremely low 77 values 
due to the method for both procedures. As a consequence the 
majority of the lines showing low 77 values from profile fitting 
should be located close to the cut-off of the b(N) distribution. 

Fig.|8]shows the measured H 1 column density versus the b- 
parameter (both logarithmic) as well as log A^h i versus the col- 
umn density ratio 77 for the combined line sample of HE 2347- 
4342 and HS 1700+6416. The right panel illustrates the corre- 
lation of 77 with the strength of the H 1 absorption. Points with 
log 77 < 1 .0 are highlighted. Their average column density is 
13.80 + 0.95. Less than 20 % of the low 77 lines have Hi col- 
umn densities < 10 13 crrT 2 , supporting the suspicion that only 
absorbers with column densities above this threshold are af- 
fected. In the N-b diagram points with log 77 < 1 .0 accumulate 
close to the lower envelope. This means, thermal broadening is 
expected to be relevant for these lines. 

We rou ghly estimate the cut-off ^-parameter by an iterative 
procedure jschave et al.1ll999l l2000t IKim et alJEool includ- 
ing only lines with crt/b < 0.5 and cr\ ogNii J log Nm < 0.5. 
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Fig. 8. Measured H i column density versus logarithmic Doppler parameter (left panel) for the combined line sample of HE 2347- 
4342 and HS 1700+6416. The solid line represents the iteratively estimated cut-off ^-parameter log b c = 0.118 • logA^i - 0.26 
(see text). Only lines with 12.5 < logA^i < 15.0 and crh/b < 0.5 and <Ti og N H J log Nm < 0.5 are considered. The right panel 
shows are the distribution of 77 with the H 1 column density. Lines with log 77 < 1 .0 are highlighted. 



According to iHui & Gnedinl Jl997l) . the tight power law rela- 
tion between the temperature and density is only valid in the 
low density IGM, which gives raise to the majority of the Lya 
forest features. Using the relation between hydr ogen column 
density and overdensity 6 given by Schave ( 2001 ), an overden- 
sity of 6 « 10 roughly corresponds to I orA^ht ~ 15.0. In order 
to avoid biases due to incompleteness jKim et alj |2002). only 
lines above the completeness limit are considered. We con- 
strain the sample used for the fit to lines with 12.5 < log iVm < 
15.0. The final sample contains about 66 % of the total lines 
and has a mean redshift of (z) = 2.58. 

For the model \ogb c = (F - 1) ■ logA^i + logfro we find 
(r - 1) = 0.118 + 0.00 5 and logfrn = - 0.26 + 0.07 in agree- 
ment with the results of Kim et al. (2002). Considering the low 
77 lines (log 77 < 1 .0), about 70 % are within 5 cr of the cut-off, 
where we adopt cr = 0.026 the standard error of the fit estimate. 
Furthermore, the mean column density ratio of all data points of 
the fit sample within log b c + 5 cr is log 77 = 1 .3 1 + 1 .08 (median 
1 .25), which is lower than the mean value of the total fit sample 
1 .56 + 1 .08 (median 1 .60). The effect gets more pronounced for 
data points located closer to the cut-off. Considering the points 
within 1 cr the median value of log 77 is 1.07. Extrapolating the 
estimated cut-off to lower and higher column densities, we find 
the same effect, i.e. the median column density ratio decreases, 
the closer the considered lines are located to the cut-off. Within 
5 cr (2cr, 1 cr) the median value is log 77 = 1.45 (1.39, 1.24). 

Fig. |2l illustrates the dependence of the inferred column 
density ratio from the distance d = |log/? c - log/>Hil to the 
cut-off in the b(N) distribution, i.e. from the contribution of 
thermal line broadening. The cumulative distribution of the me- 
dian 77 values is given for the total (logA^i < 15.0) sample as 
well as a low and high density subsample, which comprise the 
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Fig. 9. Cumulative distribution of the median log 77 with the dis- 
tance d = I log b c - log />h 1 1 from the b(N) cut-off. The me- 
dian column density ratio is estimated considering all absorbers 
within the given distance from the cut-off log b c . The solid line 
represents all absorbers with cri,/b < 0.5, (Ti os n h I log A^hi < 
0.5, and logA^i < 15.0. Additionally, the distributions of 
the low density sample containing absorbers with logA^i < 
13.0 (dotted) and the high density sample with lines 13.0 < 
logiVei < 15.0 (dashed) are indicated. 

lines with logA^i < 13.0 and 13.0 < logA^i < 15.0, respec- 
tively. It is clearly seen that the low density sample approaches 
a larger median 77 value (log 77 = 1 .93) than the high density 
sample (1.31). The low density sample reaches the plateau val- 
ues close as 98 % at d « 0.03, where only 4 % of the data points 
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are used to estimate the median. The same level, i.e. 98 % of the 
plateau value, is reached at d * 0.65 in the case of the high den- 
sity sample, where 98 % of the absorbers are included. These 
numbers illustrate that high density absorption features are sta- 
tistically more sensitive to the distance the the b(N) cut-off than 
low density lines. The cumulative distribution of the total sam- 
ple reaches a plateau value of 92 % of that of the low density 
sample, i.e. log 77 = 1.77, but follows the smoother rise of the 
high density sample. 98 % of the plateau value is reached at 

We conclude that absorption lines dominated by thermal 
broadening, i.e. located close the the thermal cut-off of the b{N) 
distribution, reveal systematically low column density ratios 77 
when applying standard analysis procedures. In consequence, 
the apparent correlation between the strength of the H i absorp- 
tion and high 77 values should be reconsidered. As illustrated 
in Fig. the column density ratio of weaker Hi absorbers 
(log N < 13.0) will be well recovered even if the assumption of 
turbulent broadening is incorrect. While weak H i lines suppos- 
edly represent an unbiased sample (represented by the low col- 
umn density sample shown in Fig.|4j, strong lines are expected 
to be biased towards lower 77 values. Hence, weak lines are ex- 
pected to show on average higher column density ratios than 
strong lines (see also Fig.|9jl- Thus, the correlation of the col- 
umn density ratio with the strength of the H i absorption can be 
explained by inadequate handling of thermal broadened lines 
by the applied analysis procedures. 

According to this result, only the low density Lya for- 
est should be used to estimate the column density ratio. 
Considering only absorbers with logA^i < 13.0 the mean 
77 value of the combined sample of HE 2347-4342 and 
HS 1700+6416 is log 77 = 1.98 + 0.81 (median 1.93). If the 
patchy zone in the spectrum of HE 2347-4342 is neglected, we 
yield log 77 = 1.89 + 0.72 (median 1.90). For the Lya forest to- 
wards HE 2347-4342, we find 1.90 + 0.74 (median 1.88), and 
towards HS 1700+6416 1.89 ± 0.71 (median 1.93). These val- 
ues are roughly 0.3 dex higher than those obtained including 
the high density absorbers. The spectrum fitting method leads 
to values roughly 0.2 dex lower than those derived from the 
low density samples but are significantly higher than the av- 
erage values for the total line samples, implying that it is less 
sensitive to the systematic errors due to thermal line widths. 
However, this is only valid on scales Az > 0.1. Variations of 
77 on smaller scales are possibly affected by methodical uncer- 
tainties due to thermal line broadening. 

Further support is provided by the photoionization models 
of the met al line systems in the spectrum of HS 1700+6416 
derived hv IFechner etafl (|2006b). The syste m at z = 2 .4965, 
for e xample is modeled using a modified Haardt & Madau 
J200ll) UV background (HM3). Doppler parameters are com- 
puted by comparing the model temperature with the observed 
^-parameter of a we ll-measured ion (for further details see 
iFechner et alJl2006bl) . For this system, the velocity ratio £ = 
^Hen/^Hi is ~ 0.50, i.e. the line width is determined by 
pure thermal broadening. The modeled column density ratio 
is log 77 = 2.35. Fitting line profiles, we assume pure turbu- 
lent broadening and find log 77 = 0.03 ±0.10 (clearly seen in 
the lower panel of Fig.|5]at z ~ 2.5). This supports the above 



argumentation that the column density ratio is severely under- 
estimated if absorbers with thermal line widths are analyzed 
assuming turbulent line broadening. Conside ring all metal line 
systems modeled by IFechner et alJ J2006tJ> . the mean veloc- 
ity ratio is £ = 0.59 ± 0.08, where the column density of the 
Hi components is on average logA^i = 15.5 ± 0.7. Thus, ac- 
cording to the photoionization models the line widths of strong 
metal line systems towards HS 1700+6416 are basically due 
to thermal broadening and can therefore explain the apparent 
correlation between low 77 values. 

5.3. Physical implications 

Even if part of the low 77 values can be explained by system- 
atical biases due to thermal broadening, part of them might 
be real. In such cases, the coincidence of low 77 values and 
metal line absorption systems, showing features of highly ion- 
ized material, suggest the vicinity of a hard radiation source 
like an AGN. Alternatively, the presence of Ovi may also 
imply the presence of a hot, collisionally ionized gas phase. 
The suspected redshifts would be z ~ 2.40, 2.48, 2.63, 2.65, 
2.69, as well as 2.82 and 2.86 towards HE 2347-4342. Further 
support for this assumption is provided by the d etected m etal 
fines of highly ionized species (see Table |2J. iBolton et alJ 
(2006) showed that fluctuations may be due to variations in 
the number, luminosity and spectral shape of the small num- 
ber of QSOs. In our case, the redshift range z ~ 2.40 towards 
HE 2347-4342 is of particular interest. It is the center of an 
apparent void in Hi, where we measure low 77 values. We find 
log 77 = 1 .07*J'H for 2.39 < z < 2.40. In the optical some weak 
absorption features with logA^i < 13.5 are detected. The red- 
shift ranges 2.376 < z < 2.384 and 2.401 < z < 2.419, where 
only Hi lines with logA^i 5= 12.7 are detected, exhibit high 77 
values (log 77 = 2.14^ ^ and 2.41+yj^, respectively). Besides 
fluctuations due to statistical properties of the ionizing sources, 
this configuration may indicate the spatial vicinity of a QSO. In 
this case the z ~ 2.40 region towards HE 2347-4342 would be 
one of th e rare cases of a classical tra nsversal proximity effect 
fe.g. lSchirber et alJ2004lEroftl2004l) . However, since there are 
no published positions of QSOs close to this line of sight for 
the required redshift range, its presence cannot be confirmed to 
date. 

Towards HS 1700+6416 the amplitude of the observed 
fluctuations is apparently smaller, in particular in the redshift 
range 2.58 < z < 2.75 where the data quality is best. This 
may be partly due to the lower noise level of the FUSE data. 
HS 1700+6416 was target of deep direct observation cam- 
paign s with the aim to i dentify galaxies close to the line o f 
sight dTeplitz et alJll998HErb et al|2003UShapley et alJ | 2005h . 
Among the galaxies identified by Shaplev et alJ (120051) there 
are three objects showing indications of quasar activity. They 
are located at z ~ 2.293, 2.333, and 2.347, which is at the low 
edge of our FUSE data. We see indeed indications of relatively 
low 77 regions in this redshift range. However, due to the poor 
data quality in this part, this coincidence should be considered 
very preliminary. 
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Fig. 10. Distribution of estimated sizes of the scales of 77 
variation towards HE 2347-4342 (left panel; z < 2.74) and 
HS 1700+6416 (right panel). The left panel also shows the dis- 
tribution of scale sizes towards HE 2347-4342 considering the 
best part of the He n only (filled histogram; 1088 < A < 1 134 A 
corresponding to 2.58 < z < 2.73). 
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Az 

Fig. 11. Scale size versus 77 value for HS 1700+6416 (circles) 
and the He n Lya forest of HE 2347-4342 (triangles). 

5.4. Spatial scales ofrj variations 

The spectrum fit method leads to estimates of the spatial scales 
of j] variations. Due to the fitting procedure the smallest de- 
tectable scales hav e a size of Az » 0. 001 which corresponds to 
the value found bv lShull et all J2004 towards HE 2347-4342. 

The resulting distributions of scales as derived by our % 2 - 
procedure (the middle panels of Figs.[3]and[5j are shown in Fig. 
[TUI Towards HE 2347-4342 only the Lya forest part z < 2.73 
of the Hen spectrum is considered, i.e. the patchy not com- 
pletely reionized part is not used. In this case scales in the range 
0.0013 < Az < 0.0473 are derived. Obviously, this sample fa- 
vors small scale variations since 69 % of the scales have a size 
of Az < 0.005, corresponding to < 4/T 1 Mpc comoving. Only 
20 % of the estimated scales are larger than Az > 0.01, corre- 
sponding to > Sh~ l Mpc comoving. Towards HS 1700+6416 
larger scales in the range 0.0020 < Az < 0.057 1 are found. For 
this line of sight 3 1 % of the scales are smaller than Az < 0.005. 

The discrepancy of the inferred scale sizes between the 
sight lines may be due cosmic variance. However, an alternative 
explanation could be the lower signal-to-noise ratio of the He n 
data of HE 2347-4342 compared to HS 1700+6416. Limiting 
the sample to the range 2.58 < z < 2.73 towards HE 2347-4342 
where the quality of the FUSE spectrum is best and the S/N is 
comparable to the HS 1700+6416 data, leads to a distribution 



similar to that towards HS 1700+6416 (indicated as the filled 
histogram in the left panel of Fig.fTol. Though the limited sam- 
ple provides only low statistics, we find that 33 % of the scales 
have a size of Az < 0.005. Large scales of Az > 0.01 are de- 
rived for 39% of the scales (46% towards HS 1700+6416). 
This indicates that part of the small scale variations might be 
mimicked by the high noise level of the He n data. In order to 
clarify this point spectra with improved S /N (at least by a factor 
of 2) would be needed. However, the results presented in this 
work indicate that a the large amplitude of the small scale fluc- 
tuations would decrease for data with a higher signal-to-noise 
ratio. 

Differences between small-scale and large-scale variations 
are illustrated in Fig.^Jwhich presents the distribution of the 77 
values with the scale size. Obviously, the scatter of the 77 values 
increases with decreasing scale size. In particular, the 77 values 
for scales Az < 0.005 towards HE 2347-4342 are extremely 
uncertain. The mean 77 value on scales Az < 0.01 for the com- 
bined sample of both lines of sight is log 77 = 2.30 ± 1.60 (me- 
dian 2.05), while scales of Az > 0.01 lead to log 77 = 1.86+0.33 
(median 1.87). This means the scatter of 77 for small scales is 
more than 1 dex larger than for larger scales. Considering both 
quasar sight lines individually leads to the same result. 

Thus, we conclude from the analysis of the spectra of 
HS 1700+6416 and the best part of the HE 2347-4342 data 
that about 30 % of the 77 variations occur on scales < 4h~ l Mpc 
comoving, while roughly 50 % of the variations have scales 
> 10 h l Mpc comoving. Large variations of 77 over several or- 
ders of magnitude apparently may occur on small scales of a 
few Mpc size, whereas there is only little scatter of less than 
+0.3 dex in 77 on large scales > 10 Mpc comoving. 

6. Summary and conclusions 

We have re-analyzed the Hen data of the QSOs HE 2347- 
4342 and HS 1700+6416, both observed with FUSE. The new 
method we applied is based on the idea to fit the observed 
high-quality optical spectrum to the noisy far-UV data glob- 
ally. This approach should minimize effects due to noise, which 
limit s the appl icability of the apparent optical depth method 
dFox etalJ [2005). and avoid the subjectivity of a profile-fitting 
procedure. Fitting the metal line cleaned Hi Lya forest to the 
He 11 data leads to log 77 = 1-79+°^ for 2.33 < z < 2.73 to- 
wards HE 2347-4342 and log 77 = 1.87+°^ for 2.29 < z < 2.75 
towards HS 1700+6416. These results are in agreement with 
the median values from the Doppler profile analysis, however, 
slightly lower than the median for lines with logA^i < 13.0, 
which we interpret as due to the neglect of thermal broaden- 
ing (see below). The mean value of the Lya forest lines with 
logA^ei < 13.0 in both quasars is log77 = 1.89 + 0.72, i.e. 
77 ~ 80. 

Strong local deviations from the average suggest that the 
column density ratio is fluctuating on smaller scales. In order 
to determine the redshift scale on which 77 fluctuates, we esti- 
mate the width of the spectral range that can be fitted with a 
minimal^- 2 . We find redshift intervals of the sizes 0.0013 < 
Az < 0.0473 in the case of HE 2347-4342. Due to the noisy 
FUSE data, the result may be biased towards small scales. For 
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the best part of the spectrum (1088 - 1 134 A), the average red- 
shift interval is Az = 0.0082 ± 0.0065 (median 0.0072, i.e. 
~ 5.9/i~'Mpc comoving). In case of HS 1700+6416, the typ- 
ical length scale is Az = 0.0140 (median) corresponding to 
~ 11 Ah" 1 Mpc comoving, but also scales up to Az = 0.0571 
are found. However, strong Hi lines, whose line widths are 
dominated by thermal line broadening, may lead to an appar- 
ent decrease of r\. Consequently, the derived scales are possi- 
bly biased towards small values. A rough estimate neglecting 
the small values leads to Az ~ 0.01 - 0.03, corresponding to 
8 - 24 /T 1 Mpc comoving. In more detail, we find that roughly 
30 % of the 77 variations take place on small scales <<\hr l Mpc 
comoving with an amplitude of about ±1.5 dex. While on large 
scales > 10 /i 1 Mpc comoving there is only little scatter of 77 
(+0.3 dex or even less). 

Towards HE 2347-4342 we detect a H i void in the range 
2.37 < z < 2.43, also showing fluctuations of 77 on scales of 
Az ~ 0.01 - 0.02 (8.5 - 17.0 A -1 Mpc comoving). The void's 
central region at z ~ 2.40 exhibits low 77 values suggesting the 
vicinity of a QSO close to the line of sight. On the other hand, 
no proximity effect can be seen towards HS 1700+6416, even 
though there are metal absorption systems showing material 
like N v and O vi supposed to be exposed to the quasar's radia- 
tion field. 

The correlatio n between H i vo ids and small 7 7 value s as 
claimed by IShull et alJ fe004 and iReimers et all J2006I) has 
been re-investigated. We find that the distribution of the column 
density ratio r\ peaks at higher values for samples consisting of 
stronger H 1 absorption lines or higher H 1 optical depth, respec- 
tively. However, simple tests with artificial data reveal that both 
analysis procedures considered here would severely underes- 
timate 77 in the case of strong Hi absorption (logA^ > 13.0) 
if the line width would be dominated by thermal broadening. 
Investigating the b(N) distribution of the line sample obtained 
from Doppler profile fitting, we find indeed that the majority 
(70 %) of lines with log 77 < 1 .0 is located within 5 cr of the 
temperature cut-off. Inversely, the average column density ra- 
tio decreases if the considered line sample is located closer to 
the cut-off. The median value is log?/ = 1.32 (1.07) for lines 
within +2cr (1 cr) from the cut-off. Photoionization mod els of 
the me tal line systems towards HS 1700+6416 (tFechner et alJ 
2006b) support this finding, since there are indications that the 
line widths of most of the modeled systems are dominated by 
thermal broadening. Thus, at least part of the low 77 values are 
due to systematic biases caused by analyzing thermal broad- 
ened absorbers under the assumption of pure turbulent broad- 
ening. Since the sensitivity of the absorbers to this effect is cor- 
related with the Hi column density, thermal line broadening 
provides an explanation of the apparent correlation between 77 
and the strength of the H 1 absorption. Consequently, only low 
density Hi absorbers provide an unbiased sample, since their 
inferred results are insensitive to the assumed line broaden- 
ing mechanism. As median column density ratios for absorbers 
with log Nh 1 < 1 3 .0, we find log 77 = 1 .93 towards both quasars 
HE 2347-4342 and HS 1700+6416. These values are about 
10 % higher in comparison to those inferred from the total sam- 
ple. 



Metal line absorption usually arises from high density H 1 
absorbers. Consequently, they lead systematically to low 77 val- 
ues. This is particularly seen in the spectrum of HE 2347-4342. 
Towards this QSO the concerned systems show highly ionized 
species like C iv and O vl The presence of high ionization lines 
suggests that the apparent correlation is not only a systematic 
effect due to the analysis method but the absorbers may be ex- 
posed to the hard radiation of local quasars. However, in or- 
der to verify this supposition, a search for QSOs in the field of 
HE 2347-4342 is required. 
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